function [ A_draw ] =restriction( phi,setup,A_draws )
%calculates a draw of the coefficient matrix when using setup.symmetry=4



end

